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Background & Motivation 


The effect aerosols have on climate and air quality is a function of their chemical composition, 
concentration and spatial distribution. These parameters are controlled by emissions, heterogeneous and 
homogeneous chemistry, where thermodynamics plays a key role, transport, which includes stratospheric- 
tropospheric exchange, and depositional sinks. In this work we demonstrate the effect of some of these 
processes on the S0 4 -NH 4 -N0 3 system using the GISS ModelE2 Global Circulation Model (GCM). 


Motivation: N0 3 aerosol is poorly constrained throughout the troposphere, especially above surface level. 
Mission: Bridge this knowledge gap with a collection of surface and airborne data and aerosol models. 
Relevant studies: Bauer at al., 2007, Bellouin et al. , 2011, Aan de Brugh et al., 2012, Hauglustaine et al., 2014, 
Paulot et al., 2015 


Objective: evaluate the GISS ModelE2 aerosol schemes and pin 
point key process either included or missing in the model 
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N 2 0 + hv- NO 


S0 2 +H 2 0-H 2 S0 3 (aq)+H 2 0 (aq) -> HSO' 3 + H + 
HSO- 3 + H 2 0 2 (aq) + H + - S0‘ 2 4 + H 2 0(aq) + 2H + 


What are the major sources of 
nitrate and sulfate precursors? 
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Figure 1: precursor sources and emissions (a compilation of references) 


How does N0 3 aerosol 
from? 

1. GAS PHASE 


h 2 o 

so 2 ->h 2 so 4 

NH, 


NO v ->HNO 


H 2 S0 4 (g) 

HNO s (g) 


2. DissolveTrfwater 

/1 

Henry'sLaw 

3. Dissociate into ions 
2H 2 0(aq) - OH(aq) + H s O + (aq) 
H 2 0(aq) - HS0 4 (aq) + H 3 0 + (aq) 
H 2 0(aq) ^ N0 3 (aq) + H 3 0 + (aq) 

NH 3 (g) f H 2 0(aq) - NH 4 + (aq) + OH(aq) 

4. Ionic solution - salt equilibrium 

HS0 4 (aq) + H 2 0(aq) • ' S0 4 2 (aq) + H s O + (aq) 
H 4 + (aq) + HS0 4 (aq) - NH 4 HS0 4 (s) 

The phase is controlled by: \ 2NH/(aq) + S0 4 2 (aq) - (NH 4 ) 2 S0 4 (s) 

1. Precursor abundance NH 4 + (aq) + N0 3 (aq) ^ NH 4 N0 3 (s) 

2. Ambient Relative humidity (Rh^ ✓ RH =water activity 

3. Temperature 

* Also in clouds 

Figure 2: The processes controlling N0 3 formation 



Setup: 

• 2° x 2.5° resolution, 40 vertical layers 

• Fully interactive trop and strat chemistry 

• Horizontal wind nudged: 6-hourly NCEP 

• SST and Ice Cover prescribed from obs 

• Emissions: 

- CMIP5 

— RCP4.5 (2005 onwards) 

— Biomass burning: GFED3 

— Agricultural NH 3 imposed seasonality 
according to solar zenith angle 


GISS Aerosol Schemes 


MATRIX : microphysics model, ◄ ►< 

tracks mixing state A/ 

OMA : bulk aerosol, includes ' 

heterogeneous uptake on dust 


EQSAM : parameterized 
thermodynamics 

ISORROPIA II : calculates 
the thermodynamics 
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particle 
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Sensitivity runs: (Result 3) we test the sensitivity of NH 3 , NH 4 , HNO 


3 > 


N0 3 to doubling and fivefold increase in agricultural NH 3 emissions. 


Data and Methods 


Monthly mean Surface data of S0 4 , NH 3 , NH 4 , NO s measured via the 
IMPROVE (USA) and EMEP-NILU (EU) networks during 2000-2010 is 
used to compare against the simulations. From the climatological 
mean (Figure 4) we adopt a regional approach (black frames in 
Figures 3,4), where the mean, standard deviation, normalized mean 
bias and correlation coefficients are calculated for the stations 
within a region along with their matching model grid boxes. 


During 2001-2011, 14 flight campaigns took place in the NH and 
measured S0 4 , NH 4 , HN0 3 , N0 3 (Figure 3). The flights were 
predominantly during spring and summer time and deployed the 
AMS instrument. With a regional approach we parse out transit 
flights and for flights within the ARC, EUSA, WUSA regions we use 
the data within the regional boundaries to calculate a campaign 
mean per model layer. We sample our simulations according to the 
flight location. 


Arctic (ARC) [55°-90°N, 60°-170°W] 

Eastern USA (EUSA) [30°-50°N, 60°-95°W] 

Western USA (WUSA) [30°-50°N, 114°-130°W] 
European Union (EU) [35°-70°N, 10°W-30°E] 



Figure 3: (above) Flight tracks of 14 flight 
campaigns used in this study 
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Results 


Is there a spatial pattern (Figure 

4)? Surface concentrations show 

high concentrations in EUSA, EU 

and low concentrations in WUSA. 

The statistics shows (Figure 5): 

* Performance is controlled by 
region more than aerosol 
scheme 

* Systematic underestimation of 
aerosols in EUSA, EU 

* Big differences for S0 4 with 
microphysics (MATRIX VS OMA) 

* Overall good performance by 
the GCM (R>0.5) 

NO ? annual cycle (Figure 6 ): 

* Seasonality is reproduced in 
EUSA, EU 

* Summer underestimations in all 


Figure 5: (below) Surface regional statistics (2000-2010); correlation coefficient 
and normalized mean bias (NMB) for the three simulations. 
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Figure 6 : 

(right) 

2000-2010 

N0 3 mean 

annual cycle. 

Error bars 

represent ± 

one standard 

deviation. 
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Figure 7: (above) 
Mean regional 
profile correspo- 
nding to ARCTAS 
spring, INTEX-A, 
CALNEX and EU- 
CCARI campaigns 
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• Model underestimation of aerosol 
concentration, especially in the 
boundary layer, more for NH 4 , N0 3 

*S0 4 : MATRIX simulations (microphysics 
included) show higher concentrations 
than the bulk scheme (OMA). 
Thermodynamic scheme makes a 
minor difference: green and red lines 
overlaying each other. 

•The model overestimates upper 
tropospheric HN0 3 : observed and 
attributed in previous studies to 
overestimation of strat-trop exchange. 

• Missing processes? (1) Heterogeneous 
chemistry sink on dust surfaces 
included in OMA, yields lower 
concentrations. (2) Homogenous 
chemistry: (i) Crustal (Mg, K, Li) and sea 
salt (Na, Cl) ions are not part of the 
thermodynamics implemented in the 
model, (ii) Thermodynamics do not 
take into account the different time 
scales associated with particle size 
distribution, (iii) Organic nitrate 
formation, especially important in the 
summer over forested areas. 

• Uncertainty in NH 3 emissions (go to 3). 
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Figure 4: (above) Mean surface concentration (2000-2010) simulated 
by MATRIX-EQSAM overlaid by measurements 
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ARCTAS spring, 

6 INTEX-A, CALNEX and 
4 EUCCARI campaigns 
plotted against NH 3 
2 emissions sensitivity 
runs. Figure 9: 
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How does increasing agricultural sensitivity runs 
4 NH ? emissions affect the system? 

•The added NH 3 acts as a sink to HN0 3 and source 
2 for NH 4 and N0 3 aerosols both at surface and above 
0 * Surface NH 3 however, appears to be already 



overestimated with normal NH 3 emissions 
Modeled and 
measured total 
NH X seem to 
match normal 
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>NH X is not partitioning properly, 
leading to underestimation of N0 3 02 
and overestimation of HNO, 
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Take away: 

> Missing aerosol mass, especially above the surface, could have important implications to aerosol radiative forcing 

> Good correlations (R>0.5) at surface in regions where seasonality is reproduced 

> HNO3 is sensitive to the heterogeneous sink - an important process to include in models 

> HN0 3 and N0 3 partitioning is strongly dependent on NH X partitioning 

> Need for more measurements: few campaigns measured NH 3 (CALNEX, TexAQ), no winter campaigns 









